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Abstract 



We prove a necessary and sufficient condition for an Abelian Sandpile Model (ASM) 
to be avalanche-finite, namely: all unstable states of the system can be brought 
back to stability in finite number of topplings. The method is also computationally 
feasible since it involves no greater than O {N^) arithmetic computations where 
is the total number of sites of the system. 
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The Abelian Sandpile Model (ASM), whose mathematical structure is first studied extensively 
by Dhar is one of the few class of models of self-organized criticality in which a lot of interesting 
physical properties can be found analytically. The model consists of a finite number of sites labeled 
by an index set /. For each site ? G /, we assign an integer hi called the local height to it. Whenever 
the local height of a site exceeds a threshold (which is fixed to for simplicity), the site is called 
unstable and it will transport some of its local heights (or sometimes called particles at that site) 
to the other sites in the coming timestep according to 

hj — > hj — Aij whenever hi > 0. (1) 

A is called the toppling matrix whose elements satisfies 

Aii>0 y le I, (2a) 

and 

A,, < V 2 ^ J. (2b) 

Sometimes, we may further require that particles cannot be created in the re-distribution process, 
which means that 

^ A,j > V z e J. (2c) 

Toppling is repeated until all sites become stable again. The whole process of toppling is collectively 
known as an avalanche. The system is driven by adding a unit amount of particles onto the sites 
randomly and uniformly after the system regains stability. 

It is, therefore, important to know the condition(s) for which the system relaxes to a stable state 
in a finite number of topplings. In fact, almost all previous studies of the Abelian Sandpile Model 
assume that such a relaxation can be taken in a finite number of steps. An Abelian Sandpile is said 
to be avalanche-finite if and only if every unstable system configuration can regain stability in a 
finite number of steps via toppling. In this paper, we prove a necessary and sufficient condition for 
avalanche-finiteness. Actually our proof works equally well on integer or real toppling matrices. 

A system configuration, stable or not, can be regarded as a point in the space where 
is the total number of sites in the system. Moreover, toppling of particles can be regarded as a 
translation in M^. If site z of a system configuration a is unstable, then the system configuration 
will become a — A, in the next timestep where Aj is the i-th row of the toppling matrix A [0,0. 

Dhar showed that the order of toppling does not affect the outcome of the final stable state in 
an avalanche. In this respect, the model is commutative. He then went on to prove that the average 
number of toppling occurs in site j given that a particle is introduced in site i, Gij, is given by 

= (A"')., (3) 

provided that A is invertible 

Now we are going to prove a necessary and sufficient condition for avalanche-finiteness. 

Proposition 1: An Abelian Sandpile Model with toppling matrix A is avalanche-finite if and only if 
there exist rzj e Z+ Vz G / such that Aj^ > for all j G /. 
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Proof: (^) Consider the unstable configuration a = (1, 1, ... , 1). Avalanche-finiteness implies the 
existence of rij G Z"*" such that a — J2i ^i^i is a stable configuration. As a result, 1 — njAj^ < 
for all J G / and hence it is proved. 

(<^=) Suppose A is not avalanche-finite, then there exists an a G which cannot regain stability 
by toppling in a finite number of steps. Let J be the set of all unstable sites at this moment. Define 

/3« = by 

V / iei 

(1) ri ifjGj 

\ if J ^ J ■ ^ ' 

Then in the next timestep, a will become a — /3A. As the toppling will never stop, we will get an 
infinite sequence of vectors by repeating the above process. Now if 1 occurs finitely many 

times in then site j will be stable eventually. Let I' be the set of all sites which become 

stable in a finite number of steps. Clearly I ^ I' or else the toppling stops in finite number of steps. 



It is easy to see that Aij = for z G /\/' and j G /'. Besides, Eq. (|2BD together with ^^^j riiAij > 
implies that ^i^ij > for all j G /\/'. But this means the local height for sites in /\/' will 

be negative after sufficiently long time, which is impossible. □ 

Proposition 2: An Abelian Sandpile Model with toppling matrix A is avalanche-finite if and only if 
det A 7^ and = (A-^).^. > for all i,j G /. 

Proof: Let Pj = {J^jj^i^j^j '■ bj > Wj} be the non-negative span of {Ajj^g/^jj}, and C = 
{Ej6/ ^i^i : &i > Vj} be the cone generated by Eq. (^) tells us that Uie/ (M+)^ = 

0. 

We consider the situation where det A 7^ first. As {Aj}j^j is linearly independent, [j-^j Pi, 
the boundary of C, disconnects M^. So either (i) C D (M+)^; or (ii) C n {R+f = 0. 

Case(i): ^ (M+)^ A'^ C (M+ U {0})^. That is, (M+)^ A'^ is an open cone in (M+ U {0})^. 
Therefore, (M"*")^ A~^ contains (infinitely many) positive integer lattice points. Hence result follows 
from Proposition 1. 

Now we consider the remaining possibility where det A = 0. As {Ajj^g/ is linearly dependent, 
one can check that Jig/ Pi = C, hence C n = 0. Thus, for any Ui G Vi G /, there exists 

j G / such that ^lUiAij < 0. By Proposition 1, the system is not avalanche-finite. □ 

With the above two propositions at hand, we are ready to prove our main theorem. 



Theorem 1: An Abelian Sandpile Model with toppling matrix A satisfying Eqs. (^a[) and (^) is 
avalanche-finite if and only if det A > and Gij = (A~^)-^- > for all i,j G /. 

Proof: By Proposition 2, it suffices to show that "the ASM is avalanche-finite implies det A > 0". 
We prove this claim by induction on N. 

Obviously, our theorem holds when = L Suppose it is true for = A; — 1. Now for N = k, 
Proposition 2 assures that det A 7^ 0. Let A' be the (A; — 1) x (A; — 1) sub-matrix formed by deleting 
the A;-th row and A;-th column of A. Using Eq. ( |2ED and Proposition 1, A' is avalanche-finite with 
k — 1 sites. So by induction hypothesis, Ckk{A) = det A' > where Ckk{A) denotes the (fc, fc)-th 
cofactor of A. By Proposition 2, (det A) ■ Ckk{A) > 0. Therefore, det A > 0. So it is proved. □ 
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Remark 1: For = 1,2,3, one can easily check by direct computation that an Abehan Sandpile 
Model is avalanche-finite if and only if det A > 0. 

Remark 2: For > 4, det A > alone is not a sufficient condition for avalanche-finiteness. Consider 
the system with A given by 



1 -2 
-1 1 







1 

-1 



'Af-4 



(5) 



with a = (1, 0, 0, . . . , 0). Direct computation verifies that the system is not avalanche-finite. 

Finally, if we restrict ourselves to the case where particles cannot be created during the toppling 
process, a simple criterion for avalanche-finiteness is found. But before we state our finding, we first 
report two important lemmas. 



Lemma 1: det A > if A satisfies Eqs. (PaD-(PcD. 
Pf: Consider the elementary row operations of 



row i\ 



row i\ 



Aa 
An 



row 1], 



(6) 



for i > 2, with resulting matrix A*^^) 
Consider 



A 



(1) 



(1) 



. Now A^^; 



(1) 



A 



(1) 

22 



A22- 



A21I lA 



12 



An 



> A 



22 



|A 



211 



for i > 2, A^^ 



A22 + A21 > 



< for i ^ j, i > 2. 



(7) 



22 



with equality holds if and only if An = IA12I and A 
Aij = A2j = for j > 2. So A^j = if and only if A^^^ 



On the other hand, if A22'' > 0, then A 
,(1) 



: IA21I. Eqs. (^) and (|2bD implies that 
0. Then det A = 0. 



(1) 



> 0, A^^^ < for i 7^ j, > 2. Also, it is easy to verify 

that A,^;^ > 0. Hence we can repeat our argument on a smaller sub-matrix until either (i) we 

get to the point where the equality holds in Eq. (|^ and hence det A = 0; or (ii) we eventually end 
up with an upper triangular matrix with positive diagonal elements, and hence det A > 0. □ 

Lemma 2: Cjj(A) > for all z, j G / if A satisfies Eqs. (|j 
cofactor of A. 



^). Here Cjj(A) denotes the (z, j)-th 



Proof: Note that Cij{A) = det H, where S is the N x N matrix with elements given by 




if JO = i and q = j 
ii p = i and q j 
if p ^ i and q = j 
otherwise 



(8) 



It is easy to check that the row sums of H are all non-negative except for the j-th row. Let k be an 



integer greater than ^ 



r&I\{j} l^i'-l' 



and consider the row operation [row j] — > [row j] + A; [row i]. 
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Then it is obvious that the row sums of this new matrix are all non-negative and hence by Lemma 1, 
Cij{A) = detE>0. □ 

Remark ^: Thus if we further require A to be invertible, then Gij = {A~^)-j > for all i,j G /. 

Corollary 1: An Abelian Sandpile Model with toppling matrix A satisfying Eqs. (^aD-(^ is 
avalanche-finite if and only if det A > 0. 

Proof: Direct application of Theorem 1 and Lemma 2. □ 

In conclusion, we found a necessary and sufficient condition for avalanche-finiteness for Abelian 
Sandpiles. In the event that particles may be created during the toppling process, avalanche- 
finiteness can be examined by looking at the determinant and the elements of the inverse of the 
toppling matrix A, a process requiring not greater than O (A^'^) arithmetic computations. Thus, our 
method is computationally practical. Furthermore, our finding is consistent with the interpretation 
that (A-i).^. is the average number of topplings occurs in site j given a particle is introduced in site 
i. 

If we further restrict ourselves to the case where particle cannot be created in toppling, our 
avalanche-finiteness testing scheme can be further simplified to a simple test in the sign of det A, 
which again can be done efficiently in the computational point of view. 
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